Three-terminal transport through a quantum dot in the Kondo regime: Conductance, 

dephasing and current-current correlations 
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We investigate tiie nonequilibrium transport properties of a tiiree-terminal quantum dot in tiie 
strongly interacting limit. At low temperatures, a Kondo resonance arises from the antiferromagnetic 
coupling between the localized electron in the quantum dot and the conduction electrons in source 
and drain leads. It is known that the local density of states is accessible through the differential 
conductance measured at the (weakly coupled) third lead. Here, we consider the multiterminal 
current-current correlations (shot noise and cross correlations measured at two different terminals). 
We discuss the dependence of the current correlations on a number of external parameters: bias 
voltage, magnetic field and magnetization of the leads. When the Kondo resonance is split by 
fixing the voltage bias between two leads, the shot noise shows a nontrivial dependence on the 
voltage applied to the third lead. We show that the cross correlations of the current are more 
sensitive than the conductance to the appearance of an external magnetic field. When the leads 
are ferromagnetic and their magnetizations point along opposite directions, we find a reduction of 
the cross correlations. Moreover, we report on the effect of dephasing in the Kondo state for a 
two-terminal geometry when the third lead plays the role of a fictitious voltage probe. 

PACS numbers: 72.15.Qm, 72.70.-fm, 73.63.Kv 
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I. INTRODUCTION 

The Kondo effect represents a distinguished example 
of strong many-body correlations in condensed matter 
physics.^ Over the last fifteen years, much effort has been 
made in understanding the implications of the Kondo ef- 
fect on the scattering properties of phase-coherent con- 
ductors. Indeed, the electric transport through a quan- 
tum dot connected to two terminals becomes highly cor- 
related when the temperature T is lowered below a char- 
acteristic energy scale given by ksTx^ At equilibrium, 
the Kondo temperature Tk depends on the parameters 
of the system, i.e., the coupling of the dot to the external 
leads due to tunneling, the dot onsite repulsion (charg- 
ing energy) and the position of the resonant level relative 
to the Fermi energy Ep- All of them can be tuned in a 
controlled wayi^ 

In a quantum dot with a sufficiently large charging en- 
ergy [U ^ kgT) and a single energy level well below Ep, 
the dynamics of the quasilocalized electron becomes al- 
most frozen. Therefore, a quantum dot can be viewed as 
an artificial realization (at the nanoscale) of a magnetic 
impurity with spin S — 1/2. At very low temperatures 
(T < Tk), charge fluctuations in the dot are suppressed 
and there arises an effective antiferromagnetic interaction 
between the electrons of the reservoir and the 5=1/2 
localized moment. Remarkably, the measured conduc- 
tance reaches the maximum value for a quantum channel 
(2e^//i) and the dot appears to be perfectly transparent 
when a small voltage eVsd is applied between the source 
and the drain contacts. Nevertheless, the coherent cor- 
related motion of the delocalized electrons forming the 
Kondo cloud can be disturbed when either the bias volt- 
age or the temperature are of the order of Tk- In such a 
case, the many-body wave function of the Kondo state is 



expected to suffer from dephasing, leading to a decrease 
in the conductance. This issue has recently attracted a 
lot of attention^i^. In this work, we mimic in a phe- 
nomenological way the effect of dephasing on the trans- 
port properties of a two-terminal quantum dot in the 
Kondo regime by introducing a fictitious voltage probe. 

Now, in the absence of dephasing, the building block of 
the Kondo resonance is a narrow peak (of width ^ ksTK) 
around Ep in the local density of states (LDOS) of 
the dot. However, full quantum-dot spectroscopy of 
the LDOS cannot be accomplished with a two-terminal 
transport setup. In particular, one cannot gain experi- 
mental access to the predicted large voltage induced split- 
ting of the LDOS when eV^d > HbTk^^^^ A way to 
circumvent this problem is by attaching a third lead, 
as demonstrated independently by Sun and Guoii and 
Lebanon and Schiller^ In subsequent laboratory work, 
De Franceschi et alm^ observed a split Kondo resonance 
by employing a slightly modified technique — one of the 
leads was replaced by a narrow wire driven out of equilib- 
rium where left and right moving carriers have different 
electrochemical potentials. 

Motivated in part by the works cited in the preceding 
paragraph, we are concerned in this paper as well with 
the nonequilibrium Kondo physics and the fiuctuations 
of the current through a quantum dot attached to three 
leads. As is well known, the investigation of the current- 
current correlations in mesoscopic conductors has been 
a fruitful area of research?^ Nevertheless, there are still 
very scarce applications to strongly correlated systems 
as the shot noise is a purely nonequilibrium property, 
and thus more difficult to treat. Hershfield^^ calculates 
the zero-frequency shot noise using perturbation theory 
in the charging energy (valid when the Kondo correla- 
tions are not large; e.g., at T > Tk)- Yamaguchi and 
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Kawamural-- choose the tunnehng part of the Hamilto- 
nian as the perturbing parameter. Ding and Ng-'^^ study 
the frequency dependence of the noise by means of the 
equation-of- motion method (also rehable for T > Tk)- 
Meir and Golub^'* perform an exhaustive study of the 
influence of bias voltage in the shot noise of a quantum 
dot in the Kondo regime. Dong and Lei^^ discuss the 
effect on the shot noise of both external magnetic fields 
and particle-hole symmetry breaking. Avishai et alm& 
calculate the Fano factor when the leads are s-wave su- 
perconductors whereas the case of p-wave superconduc- 
tivity is treated by Aono et alm^ The authors^ examine 
the behavior of the Fano factor at zero temperature when 
the formation of the Kondo resonance competes with the 
presence of ferromagnetic leads and spin-flip processes. 
Lopez et al^ make use of the two-impurity Anderson 
Hamiltonian to address the shot noise in double quan- 
tum dot systems. To the best of our knowledge, a study 
of the current fluctuations in a multiprobe Kondo impu- 
rity is still missing. This is the gap we want to fill here. 

In mesoscopic conductors, Biittiker^ shows that the 
sign of the current cross correlations depends on the 
statistics of the carriers. It is positive (negative) for 
bosons (fermions) due to statistical bunching (antibunch- 
ing). This statement is based on a series of assump- 
tions (e.g., zero-impedance external circuits, spin inde- 
pendent transport, normal thermal leads). Positive cor- 
relations can be found if these conditions are not met 
(see Ref. for references on this subject). Here, we 
just mention a few studies based on structures involv- 
ing a quantum dot. Bagret and NazaroAiS& consider a 
Coulomb-blockaded quantum dot attached to paramag- 
netic leads whereas the ferromagnetic case and the spin- 
blockade case are treated by Cottet et alm^ Borlin et 
alm^ and Samuelsson and Biittiker^® examine the cross 
correlations of a chaotic dot in the presence of a super- 
conducting lead. In the spin dependent case, Sanchez et 
alm^ find that the sign of the cross correlations is affected 
by Andreev cross reflections. In the context of quantum 
computation, measuring current cross correlations have 
been shownSi to yield a indirect identification of the ex- 
istence of streams of entangled particles. Therefore, the 
cross correlations are a valuable tool in characterizing the 
electron transport in phase-coherent conductors. 

In this work, we consider electron transport through 
a strongly interacting quantum dot attached to three 
leads. Section |n] explains the theoretical framework 
(slave-boson mean-field theory) we use to compute the 
conductance and the current-current correlations. We 
show that the expressions for the cross correlations may 
be inferred from scattering theory applied to a Breit- 
Wigncr resonance with renormalized parameters. Sec- 
tion ^O] is devoted to the results. First, we assume that 
the third lead is a fictitious voltage probe and investi- 
gate the effect of dephasing with increasing coupling to 
the probe. Then, we consider that lead as a real probe 
and relate the differential conductance measured at one 
electrode with the LDOS of the artificial Kondo impu- 



rity. We show next that the sign of the cross correlations 
of the current is negative, as expected from the fermionic 
character of the Kondo correlations at very low temper- 
ature. Moreover, we discuss the effect of bias voltage, 
external magnetic fields, and spin-polarized tunneling in 
the cross correlations. We finish this section with an in- 
vestigation of the effect of spin polarized transport in the 
shot noise. Finally, Sec. IIVI contains our conclusions. 



II. MODEL 

We model the electric transport through the quantum 
dot using the Anderson Hamiltonian in the limit of large 
onsite Coulomb interaction U oo. This way we ne- 
glect double occupancy in the dot and the Hamiltonian 
is written in terms of the slave-boson languages^ 



+ E(^^"4a.&V.+H.C.) 

ka.a 

+A(6t6 + ^/t /,-!), 



(1) 



where c^,^^ {ckaa) is the creation (annihilation) operator 
describing an electronic state k with spin a — {T, i} and 
energy dispersion Ekaa- in the lead a — {1,2,3}, £o<t is 
the (spin-dependent) energy level in the dot and Vka is 
the coupling matrix element. The original dot second- 
quantization operators have been replaced in Eq. Q by 
a combination of the pseudofermion operator /o- and the 
boson field b. Hopping off the dot is described by the 
process c\^^b^ /„: whenever an electron is annihilated by 
fa, an empty state in the dot is created by 6^ and then 
'^laa generates an electron with spin a in the conduction 
band of contact a. The boson operator b (6^) may be re- 
garded as a projection operator onto the vacuum (empty) 
state of the quantum dot. To make sure that a state with 
double occupancy is never created, a constraint with La- 
grange multiplier A is added to the Hamiltonian. 

The current operator that yields the electronic fiow 
from lead a is given by 



(2) 



where Na = J2kcr 'Aaa'^kaa- The general form of the 
power spectrum of the current fiuctuations readsS 



({/„(r),/MO)})-(/„)(/^) , (3) 



dr e" 



5Ia — la ^ la describing the fluctuations of the current 
away from its average value la — (la)- We are inter- 
ested in the zero- frequency limit of Saiii^^j)- Since the 
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FIG. 1: The system under consideration. The central island 
is a resonant level coupled to three leads. The level may be 
shifted through a capacitative coupling to a gate. In the limit 
of a vanishingly small capacitance, double occupancy in the 
dot is forbidden and Kondo effect can arise. The current- 
current cross correlations are measured between leads 2 and 
3. 



energy scale ksTx in typical experiments is of the or- 
der of 100 mK, the frequencies should be w < 2.4 GHz. 
Moreover, we shall work at T = (see below) so that the 
current will fluctuate due to quantum fluctuations only 
(we disregard thermal fluctuations). 



A. Mean- field approximation 

The mean-field solution of the Hamiltonian consists 
of considering the effect of the boson in an averaged way, 
replacing the true operator h{t) by its expectation value 
{b(t)). Within this approximation the Hamiltonian de- 
scribes nonintcracting quasiparticles with renormalized 
couphngs: Vka\/\b\ ~* ^ka- The theory is then suitable 
for studying the Fermi-liquid fixed point of the Kondo 
problem (i.e., at T <C Tk) in which the averaged occu- 
pation in the dot is always 1 . The dominant fluctuations 
in the system are those associated to spin. 

The stationary state of the boson field is determined 
from the t oo limit of its equation of motion using the 
Keldysh technique for systems out of cquilibrium:'^^'^^ 



(4) 



kaa 



where Gf^ ,^^^{t,t) = j(cL^(0/o-(0) is the lead-dot 
lesser Green function. Next, we take into account the 
constraint: 



(5) 



G^ J {t,t) = i{fl{t)fa{t)) being the dot lesser Green 
function. It yields the nonequilibrium distribution func- 
tion in the dot. 



In evaluating the above Green functions we need the 
coupling strength given by Ta.a{e) = T^Lfe l^feaP<J(e - 
Ekaa)- In the wide band limit, one neglects the energy 
dependence of T and the hybridization width is taken as 
Taa — ^aa^Ep) for —D < E < D (D IS thc high-eucrgy 
cutoff). We notice that in thc presence of Kondo cor- 
relations thc lifetime broadening becomes renormalized 
TatT fatr = ^J2k \^ka\'^^i'^ " Sfcacr) and the bare level 
Eocr is shifted to eocr = £ocr + A. We can now give the 
full expression of the Fourier-transformed lesser Green 
function: 



(e) = 2z 



X/Q ^aafa (e) 

(e-£o.)2 + rr 



(6) 



where T^- — is the total hybridization width 

per spin and /a(e) = OifJ-a ~ e) is the Fermi function at 
zero temperature of lead a with electrochemical potential 
/ia = Ep + eVa- On the other hand, G^^ kaai'^) '-^^ ^® 
cast in terms of G^^ (w) with the help of the equation 
of motion of the operators and then applying the ana- 
lytical continuation rules in a complex time contour^ 
Therefore, we obtain a closed system of two nonlinear 
equations [Eqs. Q and ©J with unknowns |6p and A to 
be found self-consistently. 

From the precedent arguments and Eq. |(5J we can eas- 
ily establish an expression for the expectation value of the 
electric current: 



I3a 



def^^{e)[Ue)-Me)], (7) 



which has exactly the same transparent form as the 
Landauer-Biittiker formula^ in the two channel (one per 
spin) case applied to a double-barrier resonant-tunneling 
system: 



(e-£o.)2+f2 



(8) 



which has a simple Breit-Wigner lineshape. For the same 
reason, the quasiparticle density of states is a Lorentzian 
function centered around the Fermi level (the Abrikosov- 
Suhl resonance). This result is expected since we are 
dealing with a Fermi liquid but we stress that the physics 
it contains should not be confused with a noninteracting 
quantum dot since: 

(i) T depends implicitly on |6p and A, and it must then 
be self-consistently calculated for each set of parameters: 
contact voltages {V^}, magnetic field Az — £oT ~ ^oi, 
gate voltage So(Vg), and lead magnetization. 

(ii) T is renormalized by Kondo correlations (as the 
bare F and eo are), 

(iii) T has a nontrivial dependence on the bias voltage. 
All these features give rise to a number of effects 

that are not seen in a noninteracting resonant-tunneling 
diode. There are many instances: regions of negative 
differential conductance in the current-voltage charac- 
teristics of a double quantum dot^ a crossover from 
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Kondo physics to an antiferromagnetic singlet in the two- 
impurity problem an anomalous sign of the zero-bias 
magnetoresistancef^etc. Below, we shall discuss another 
example without counterpart in a noninteracting Breit- 
Wigner resonance: When the Kondo peak splits due to 
a large bias voltage. 

B. Current-current correlations 

We consider now the current fluctuations given by 
Eq. lO at zero frequency S'q^(O). To simplify the no- 
tation we introduce G'o(w) = Gf^,f^{uj) as the dot Green 
function. After lengthy algebra, we have 

4e2 /■ - ^ ^ ^ ^ 

SapiO) = I deVaVplG^G^ - GqGq fa 

+6*0 Go(1-//3)-G^Gq(1-/q)+G'oG5'//3-GqGo/q(1-/;3) 
- G5G5/;3(1 - /.) - *%(G<(1 - - G>/„)] . (9) 

This formula (or variations of it) has been already em- 
ployed in the literature. Wei et al?^ prove it using the 
Fisher-Lee-Baranger-Stone relation^S to write the scat- 
tering matrix elements in terms of the retarded Green 
function of the dot, Gq. Dong and LeiiS and Lopez et 
alMiP find it in Kondo problems within a slave-boson 
mean-field framework. Actually, in Ref. 23, it is shown 
that the shot noise in a two-terminal geometry reads 
S ~ T(l — r), i.e., the well known result for the par- 
tition noise but with renormalized transmissions. Souza 
et al^ calculate the noise of an ultrasmall magnetic tun- 
nel junction by means of Eq. (jSJ within a Hartrec-Fock 
framework. In general, we can say that Eq. ® IS consis- 
tent within mean-field theories. However, some caution 
is needed if one wishes to go beyond a mean-field level. 
In deriving Eq. ©, one needs to apply Wick theorem, 
which is valid only for noninteracting (quasi)-particles. 
More specifically, one finds terms that read: 

(cL.w/.WcU'(o)/.'(o)) = 

+ (4a.W/<.'(0))(/<xWcI/3.'(0))- (10) 

The first term in the left-hand side corresponds to dis- 
connected diagrams that cancel out the term {I a) {I 13) of 
Eq. ||2Jl whereas the second term contributes to Eq. I^. 
Therefore, the particular Hamiltonian has to be cast first 
in a quadratic form. Zhu and Balatsky'*^ incorrectly state 
that Eq. (O takes into account the many-body effects. 
Also, it is not clear how this formula is inferred within the 
equation-of-motion method employed by Lii and Liu4^ 
In our case, the mean-field approximation is known to 
be the leading term in a 1 /N expansion;^ where N = 2 
is the spin degeneracy. Therefore, we neglect the fluctu- 
ations of both the boson fleld {5h = 0) and the renormal- 
ization of the resonant level {5\ = 0),iM2, which could 



be calculated in the next order. This is valid as long as 
we restrict ourselves to the Fermi-liquid fixed point of the 
Kondo problem. We are not aware of real l/N correction 
calculations of shot noise. Although Meir and GolubA^ 
perform a noncrossing approximation (NCA), they just 
substitute the NCA propagators into Eq. 0, with the 
limitations exposed above. 

The current-current correlations can be deduced ei- 
ther using Eq. ^ or using the scattering approach for 
the multiterminal case (see Ref. |2^ . The latter formal- 
ism amounts to replacing the bare parameters by the 
renormalized ones^i. We consider the illustrative case of 
having different electrochemical potentials in two leads, 
^J'a 7^ M/5 (e.g., a — 2 and /3 = 3) at zero temperature. 
We find 

2g2 r 

523(0) = — / deTv{sl^S2SslsS^^){f^-fa){f5-h) , 

7,5 

(11) 

where Sap is the renormalized scattering amplitude of a 
Breit-Wigner resonance: 

= 5o.p - \ , - . (12) 

In Eq. (|11|) the trace Tr(...) is over spin indices. The 
Fermi functions fa and are arbitraryr^** Choosing fa = 
fb = fa, we obtain 

^23(0) = ^ / m2T^%[h - fa? 

+ ^22^3^2 [/2 - /3]' + 2f f2f2-3 [/i - /g] [/2 - /3] } , (13) 

where i?22 is the refiection probability (in general Raa = 
1 — X^/j^a^a)- Notice that generally one cannot write 
the multilead current -current correlations only in terms 
of transmission probabilities as in Eq. (|13l) . This was 
firstly pointed out by Biittiker,^^ suggesting the appear- 
ance of exchange effects in noise measurements. Here, 
since we are dealing with a (renormalized) Breit-Wigner 
resonance, exchange corrections due to phase differences 
do not play any role. 

For completeness, we give now the formula for the shot 
noise, i.e., the current-current correlations measured at 
the same lead (e.g., lead 1). Following the way of reason- 
ing that led to Eq. (O we obtain 

+ f^^Rl, [A - /3]2 + n^T^Ah - f^n ■ (14) 

III. RESULTS 

In the following, we present results obtained by self- 
consistently solving Eqs. I@J and lO for each bias volt- 
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age. The rest of parameters is changed in the next sub- 
sections. Throughout this work, we have checked that 
current conservation (/i + /2 + -^3 = 0) is fulfiUed.^^ 

Tunnehng effects are incorporated at aU orders since 
at equihbrium the Kondo temperature is found to be 



fcsT^. = f = i?cxp(-7r|£o|/2r), 



(15) 



which is clearly a nonperturbative result. In Eq. H15() 
r = X]a=i is the total hybridization broadening. The 
reference energy will be always set at Ep — and the 
energy cutoff is D = lOOF. The bare level is Eq = — SF, 
deep below Ep to ensure a pure Kondo regime. 



A. Dephasing 

Before turning to the determination of current cross 
correlators, we briefly discuss with an application the ca- 
pabilities of t/iree-terminal setups to illustrate some dif- 
ficult aspects of the physics of the two-teiminal Kondo 
effect. As mentioned in the Introduction, we investigate 
the action of a fictitious voltage probed (say, lead 3) in 
order to simulate decoherence effects on the formation 
of the Kondo resonance between leads 1 and 241 These 
contacts play the role of source and drain, respectively. 
The voltage probe model^^ describes decoherence since 
an electron that is absorbed into the probe looses its co- 
herence. The exiting electron is replaced by an electron 
(with an unrelated phase) injected by the probe. 

At low temperatures the principal source of dephas- 
ing is due to quasi-elastic scatteringi^ We consider then 
a voltage probe that preserves energyi^ The current 
through the voltage probe is zero at every energy e. Thus, 
from Eq. |(7J) the distribution function at the probe reads 



. . . ri3(e)/l(e)+T23(6)/2(6) 



(16) 



We have to insert this result into Eqs. and ||SJ) and 
solve self-consistently for the hybridization couplings F 
and the resonance level sq in the presence of quasi-elastic 
scattering for each value of the applied bias voltage. Then 
we compute numerically the differential conductance G = 
dl/dVsd, where / = /i = —I2 and Vsd = V\ — Vi- 

Figurel^a) shows G for different values of the coupling 
to the probe (we set F2 = Fi) . For Fs = we obtain 
the well known zero-bias anomaly, which arises from the 
formation of the Kondo resonance at 14d = 0. As F3 
increases we observe a quenching of the Kondo peak. The 
degree of the conductance suppression depends on the 
coupling to the probe. At each bias, /is (which has to 
be self-consistently calculated) adjusts itself to fulfill the 
condition of zero net current at each energy e. Hence, F3 
is a phenomenological parameter that includes dephasing 
processes present in the quantum dot. To see this, we 
can write down the current through, say, lead 1, using 




FIG. 2: (a) Differential conductance Gn versus bias voltage 
Vi as a function of the bare coupling Vz to the voltage probe 
(reservoir 3) for Fi = Y-2 and eo = —OF. (b) Linear conduc- 
tance 011(0) showing the reduction of the peak in (a) versus 
the coupling to the voltage probe. The dots are numerical 
results where the line corresponds to an analytical formula 
(see text). 



Eqs. and lini): 



e^FiF2 



^Fi 



deAo(e)[/i(e)-/2(e)], 



(17) 



where A^ie) = — ImGQ(e)/7r is the LDOS in the dot. 
Equation (|17|l has the form of a formula for a two- 
terminal current with G'o(e) = [e — £0 + «(?! + ^2 + 
Fs)]"^. It is straightforward to show that a nonzero F3 
leads to deviations of Eq. (|17|l from the unitary limit. 

In Fig. Efb) we plot the linear conductance G = 
G{Vsd = 0) as a function of F3/F1 from the results found 
numerically. At zero bias we can find from Eq. 117|) an 
analytical expression for the reduction of the peak: 



2e2 



h 2-HF3/F1 



(18) 



It is shown in Fig. Olb) (full line). In the limit of 
Ta/ri ^ 1 a similar expression for the reduction of the 
peak was found by Kaminski et al.,'^ the source of deco- 
herence being an ac voltage applied to the dot level. 



B. Multiterminal conductance 

From now on, we consider lead 3 as a real electrode 
with tunable voltage V3. We set V3 = V2 = and vary 
the tunneling coupling F3. The self-consistent results of 
Eqs. (@J) and Q are inserted in Eqs. {Tj) to calculate the 
differential conductance through lead 1: Gn = dli/dVi. 
FigureOl^a) shows Gn as a function of Vi. At F3 = the 
conductance at Vi = achieves the unitary limit as in 
the two-terminal case. With increasing the coupling to 
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FIG. 3: (a) Differential conductance Gu versus bias voltage 
Vi for eo = -6r (T]^ = 8 x lO^^F). (b) Dependence of the 
Kondo temperature on Vi. 



third lead, 6*11(0) decreases. For equal tunnel couplings 
(Fi = Fa = Fa = r/3), 6*11(0) does not reach 1 (in units 
of 2e^ / h) but instead Gii(O) = 8/9, in agreement with 
Ref. I51L This is an immediate consequence of having 
three leads with identical couplings. Interestingly, the 
Kondo temperature of Fig. |SIb) does not vanish abruptly 
for Vi = 2T]^, as known in the two-terminal case (see the 
case F3 = 0). This is an important result as it implies 
that Kondo correlations survive at large voltages. The 
effect is reminiscent of the situation found by Aguado 
and LangrethSi in tunnel-coupled double quantum dots, 
though the physical origin is clearly distinct. 



C. Sign of the current cross-correlations. 
Comparison with a noninteracting quantum dot 

We now focus on the current-current correlations of 
the current for V3 = V2 — and equal couplings Fi = 
F2 = F3 = F/3. Later, we shall allow for nonzero voltage 
differences between leads 2 and 3. In Fig.^a), we show 
the cross correlator 5*23(0) obtained from Eq. H13|) . As 
expected, 5*23 is zero for Vi = and negative elsewhere. 
This reflects the fermionic nature of the quasiparticles. 
For comparison, we plot in Fig. 01b) the corresponding 
5*23 for a noninteracting resonant double-barrier struc- 
ture with the level at Ep (of course, for Eq = — 6F the 
spectrum 5*23 is always very small as the transmission 
is). In this case, the physics is governed by the bare cou- 
pling F.^'^ On the contrary, in the Kondo problem the 
dominating energy scale is Tr-. Qualitatively, Fig. ^a) 
and^Jb) look the same until Vi ^ 2Tk- The cross cor- 
relator in the Kondo case increases with voltage while in 
the noninteracting case 5*23 saturates at large voltages. 
It is easy to show that the saturation value is given by 
-87r/81 ~ -0.31 (in units of Ae'^T/h). The reason for 
the increase of 6*23(0) in Fig. 2Ja) is that Tk is voltage 
dependent unlike the bare F, even in the wide-band limit. 
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FIG. 4: (a) Current-current cross correlation measured in 
leads 2 and 3, ^23(0), as a function of the bias voltage in 
the injecting lead, Vi . Kondo correlations involve an increase 
of 5*23(0) for voltages larger than 2Tk- (b) Same as (a) for 
a noninteracting quantum dot with a resonant level exactly 
at Ef- (c) and (d) correspond to the Fano factor 723 as 
a function of voltage for the interacting and noninteracting 
case respectively. 



In particular, the current-voltage characteristics shows a 
region of negative differential conductance in the Kondo 
case [see Fig.|3Ia)] whereas it reaches a constant value at 
large voltages for an noninteracting quantum dot. 

To avoid effects due to moderate biases, in what follows 
we shall concentrate on a normalized 6*23. We define the 
Fano factor of 6*23 as 



723 



'23 



2e^7^ 



(19) 



If the scattering region were a simple barrier of transmis- 
sion T, 723 would be simply —1. This number changes 
when the system under consideration is a quantum dot. 
In Figs. Ela) and (b), we plot 6*23 for the Kondo and 
the noninteracting case, respectively. Their correspond- 
ing Fano factors are shown in Figs.^Jc) and (d). We see 
that 723 has a minimum at Vi = 0. Analytically, we find 
723(0) = —4/9 ~ —0.44, which is in excellent agreement 
with the numerical result. Likewise, we can assess the 
limit of 723 at very high voltages [Vi ^ T^). We get 
723 = —2/9 ~ —0.22. As observed, both curves tend to 
this value, though for a noninteracting quantum dot it is 
more quickly due to the independence of F on the bias 
voltage. 



D. Effect of the nonequilibrium splitting on the 
current-current correlations 



Now we turn to an exciting case. Consider the follow- 
ing bias configuration: V2 = — V3 ^ and determine the 
differential conductance Gu as a function of Vi . The case 
V2 = — V3 = has been treated before. However, due to 
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FIG. 5: Differential conductance Gu versus bias voltage Vi 
for different values of the voltage difference AV = IV2 — V3I > 
2Tt. 



the fact that the boson field never vanishes, we can study 
the situation AV = IV2 - V3I > 2T^. As remarked in the 
Introduction, it has been argued^^'^^ and experimentally 
observedi^ that in a three-lead geometry the sphtting of 
the Kondo resonance due to voltage is visible , unlike the 
two-terminal case. Moreover, in Refs. |l2l5l| it has been 
noticed that the conductance Gu is not sensitive to the 
strength of the coupling to the third lead, showing al- 
ways a two-peak structure. Of course, only when the 
third lead is weakly coupled to the dot Gu is a measure 
of the LDOS. But since we are interested in the trans- 
port properties of the system, our choice of equal coupling 
constants does not affect the results for the conductances 
and the current-current correlations. 

In Fig. |31 we plot the behavior of the differential con- 
ductance Gil. At Ay = we obtain the zero-bias 
anomaly of Fig. Ola). As AV increases, Gu is split 
at V2 ~ T^. Both splitting peaks lie at Vi ^ V3 and 
Vi V2, i.e., when a pair of electrochemical potentials 
are aligned. It is also at those points where the Kondo 
temperature is larger. We emphasize that this effect has 
no similitude in the electronic transport through a non- 
interacting quantum dot. Still, a mean-field theory of 
the Kondo effect as presented here is able to capture 
this physics. At the same time that the splitting in Gu 
develops, the height of the peaks decreases, suppressing 
the zero-bias anomaly, although not so strongly as in the 
experimenlii^ due to the absence of inelastic scattering in 
this case. 

We now use Eq. H13() to calculate the cross correla- 
tions between leads 2 and 3. The results are presented 
in Fig.inja). The dependence of S23 on voltage is rather 
asymmetric, hindering the observation of a clear indica- 
tion due to the voltage induced splitting. The asymme- 
try is caused by the third term of the right-hand side 
of Eq. (|13|) , which is not symmetric under the operation 
Vi ~~V\ when AV > 0. That is the reason why we 
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FIG. 6: (a) Cross correlations of the current measured be- 
tween leads 2 and 3 for the case treated in Fig. |K| (b) Same 
as (a) for the shot noise in lead 1. 



next consider the shot noise in lead 1 6*11, which is an 
even function of the applied V\ . 

In Fig.EIb) we plot the results of Eq. H14|) . We observe 
that 5ii at Vi = is nonvanishing with increasing Al^, 
causing a divergence of the Fano factor. This is not re- 
lated to the Kondo physics but with the fact that the lead 
1 at Vi = acts as a voltage probe with zero impedance 
since the net current flowing through it is zero. Including 
the fluctuations of the potentials would probably cancel 
out the divergence. A consequence of Kondo physics is 
that the minimum at Vi = turns into a maximum. This 
occurs when the splitting in Gu is sharply formed [see 
FigE]. 



E. Spin dependent transport and current cross 
correlations 

So far we have assumed spin- independent transport. 
Let us go back to the bias configuration of Sees. IIII Bl 
and nil 01 iVo = VJj = 0) and focus on the spin-dependent 
transport properties. It is customary in the theoretical 
studies of spintronic transport to take into account the 
influence of external magnetic flelds and ferromagnetic 
electrodes, among other parameters."''^ Firstly, we shall 
change the external Zeeman field and then enable the 
presence of spin-polarized tunneling. 



1. Magnetic field 

We assume that the leads are paramagnetic and that 
the magnetic field is applied only to the dot, resulting in 
a Zeeman gap of the bare resonant level: Az = Eot ~ 
eoi- It is well known that, as a consequence, the Kondo 
resonance is split when Az ~ T^i 

Figure [7Ja) shows the differential conductance Gu for 
different values of the Zeeman field. The conductance is 
split and quenched with increasing A^, as expected. In 
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FIG. 7: (a) Differential conductance Gn versus Vn as a func- 
tion of the Zeeman term Az for V2 = Vi = 0. (b) Same as 
(a) for the Fano factor of the cross correlator, 723. 



Fig. mb), we depict the Fano factor of the cross correla- 
tor 723. It exhibits a very interesting feature. Due to the 
splitting of the Kondo peak, the minimum of the cross 
correlator at Vi = becomes a local maximum, result- 
ing from the suppression of the Kondo effect. However, 
this change occurs before the splitting of the conductance 
Gil. Therefore, measuring the shot noise provides new 
information in this case. The presence of the splitting 
would be detected in an experiment more precisely by 
means of the shot noise. The underlying reason is that 
the form of Eq. H13() differs from that of the current which 
is basically proportional to T12 alone, see Eq. 0. As a 
result, the width of the Gn resonance is a bit larger than 
the 723 antiresonance and the former is then more robust 
than the latter against the application of magnetic fields. 



2. Ferromagnetic leads 

There has recently been considerable debate about the 
influence of ferromagnetic leads in the Kondo physics of 
a quantum doti ^'^1^^1^^1^'' In the preceding subsection, it 
was clear that an external magnetic field alters the real 
part of the quantum-dot self-energy, breaking the spin 
degeneracy. In the case of spin polarized tunneling, the 
situation is more subtle.??^ When the magnetic moments 
of the contacts are aligned along the same direction, the 
density of states of the localized electron undergoes a 
splitting if particle-hole symmetry is broken^ Recent 
transport experiments with Cgo molecules and carbon 
nanotubes have addressed this regime. ^^■^'^ However, in 
our case the dot is in the strong coupling limit and the 
Kondo effect is pure in the sense that no charge fluctu- 
ations are allowed. Thus, no splitting is expected in the 
differential conductance. 

In Fig. ISJa), we show the cross correlator 723 for dif- 
ferent values of the lead magnetization in the parallel 
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FIG. 8: Fano factor of the cross correlator, 723 vs Vi/T^ for 
different lead magnetization when V2 = V3 = 0. (a) Parallel 
alignment between the magnetizations of the leads with spin 
polarizations: pi = p2 = P3 = p . (b) Antiparallel case with 
Pi = — P2 = — P3 = p. 



case. This means that pi = P2 = = P, where pa is 
the spin polarization of lead a. Ferromagnetism in the 
leads arises through spin-dependent densities of states 
t'cKT(e) = J2k ~ ^kaa)- Hcuce, the linewidths become 
spin dependent: T^a = (1 '^Pa)^on where +(-) corre- 
sponds to up (down) spins. We prefer to restrict pa to 
small values as strong magnetizations would require a 
proper treatment of the reduction of the bandwidth D. 
We observe that 723 is rather insensitive to changes in 
p in the same fashion as Gn is in the Fermi-liquid fixed 
point^^. Only at moderate polarizations {p — 0.6) we see 
that the dip in 723 gets narrower because the Kondo tem- 
perature decreases as p increases^SiS. In addition, 723 is 
always negative in contrast to the results obtained in the 
Coulomb blockade regime where 723 can take positive 
values^!. When the spin-flip scattering rate is smaller 
than the tunneling rate, 723 can be positive. However in 
the Kondo regime this condition is never met since the 
rate of spin flip scattering ~ 1/Tk is always much longer 
than the tunneling rate ^ l/F. Figure|SIb) is devoted to 
the antiparallel case: pi — — p2 — ~P3 — P- Accordingly, 
723 is lifted with increasing lead polarization since the 
conductance peak decreases with increasing p (roughly. 



with a factor 1 



p- 



2\22 



IV. CONCLUSION 

In summary, we have investigated the Kondo temper- 
ature, the differential conductance and cross correlations 
of the current when three leads are coupled to an arti- 
ficial Kondo impurity in the Fermi-liquid fixed point of 
the infinite-?/ Anderson Hamiltonian [T <^ Tk)- We 
have performed a systematic study of the properties of 
the cross correlators when dc bias, Zeeman splittings, and 
ferromagnetic leads influence the nonequilibrium trans- 
port through the quantum dot. Our most relevant result 
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is the behavior of the shot noise when there arises a volt- 
age induced spUtting in the quantum dot. 

In addition, we have studied the current of a two- 
terminal quantum dot attached to a voltage probe. We 
have shown that increasing the coupling with the probe 
induces a quenching of the Kondo peak. Despite the sim- 
plicity of this approach, it gives rise to results that are 
in agreement with more sophisticated models,i^iSi though 
the precise processes responsible for the decoherence need 
still to be derived from a microscopic model. 

We have not exhausted all the possibilities that the 
model offers and more complicated geometries with 
appealing results can be envisaged. One could ad- 
dress the situation with two injecting and two receiving 
leads, which could give rise to Hanbury Brown- Twiss-like 
effects. We expect that phase related exchange terms 
will arise especially at higher temperatures (T > Tk), 
when the singlet state between the localized spin and the 
conduction electrons is not yet well formed. We believe 
that in the presence of spin-polarized couplings due to 



ferromagnetic leads, bunching effects will be enhancediS^ 

Improvements of the model should go in the direction 
of including fluctuations of the boson field and of the 
renormalized level. However, we do not expect large de- 
viations from the results reported here when T <C Tk- 
These fluctuations will evidently become important as 
temperature approaches Tk- Experimentally, our predic- 
tions can be tested with present technology such as GaAs 
quantum dots^'^ or carbon-nanotube nanostructuresiS4 
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